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' ABSTRACT. A new derivation of the propagation direction of light is given for a 3-parameter family 

^ ^ , of static, spherically symmetric space-times within the post-post- Newtonian framework. The emitter and 

' the observer are both located at a finite distance. The case of a ray emitted at infinity is also treated. 

o' 



1. INTRODUCTION 



^ ' The aim of this work is to present a new calculation of the propagation direction of light rays in a 3- 

parameter family of static, spherically symmetric space-times within the post-post-Newtonian framework. 
^ • Rather than deriving the results from an integration of the geodesic equations, we obtain the desired 

bi), expressions by a straightforward differentiation of the time delay function (see, e. g., Teyssandier & 

Le Poncin-Lafitte 2008 and Refs. therein). This study is motivated by the fact that any in-depth 
i discussion of the highest accuracy tests of gravitational theories requires to evaluate the corrections of 

^ ' order higher than one in powers of the Schwarzschild radius (see, e.g., Ashby & Bertotti 2010 for the 

04 , Cassini experiment). Even for the Gaia mission, a discrepancy between the analytical post- Newtonian 

solution and a computational estimate has recently necessitated a thorough analysis of the post-post- 
Newtonian propagation of light (see Klioner & Zschocke 2010 and Refs. therein). 



2. LIGHT DIRECTION IN SPHERICALLY SYMMETRIC SPACE-TIMES 

The gravitational field is assumed to be generated by an isolated spherically symmetric body of mass 
M. Setting m — GM /c? , the metric is supposed to be of the form 



(1) 



where r = y/Sijx'^x^ , (3 and 7 are the usual post-Newtonian parameters, and e is a post-post-Newtonian 
parameter (/3 = 7 = e = 1 in general relativity). We put x*^ = ct and x = (cc*), with i = 1,2, 3. 

Consider a photon emitted at a point Xa at an instant and received at a point Xb at an instant 
Ib ■ The propagation direction of this photon at any point x of its path is characterized by the triple 

L={h/lo) = {h/lo,l2/lo,l3/lo), (2) 

where Iq and h are the covariant components of the vector tangent to the ray, i.e. the quantities defined 
by la = gapdx^ / d\, gap denoting the metric components and A an arbitrary parameter along the ray. 

Denote by 1^ ^^^d 1-b the expressions of Z at points Xa and Xb, respectively. In any stationary space- 
time, these triples can be derived from the relations (see Le Poncin-Lafitte et al. 2004) 

li\ _ dT{xA,XB) f h\ _ _ dT{xA,XB) , . 

IoJa~' ' Vo^~ ' dx^ ' 

where T{xa,Xb) is the expression giving the travel time of a photon as a function of Xa and x^: 

ts - tA = T{Xa,Xb). (4) 



For the metric (HJ, T{xa,Xb) is given by (see, e.g., Teyssandier & Le Poncin-Lafitte 2008): 
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Substituting for T{xa,Xb) from Eq. ^ into Eqs. yields and 1^ as Hnear combinations of 
and rig. However, it is more convenient to introduce the unit vector Nab defined by 



Na, 



(7) 



and the unit vector Pab orthogonal to Nab defined as OH /\OH\, H being the orthogonal projection of 
the center O of the mass M on the straight line passing through Xa and Xb, that is 



= AT, 



Using Eqs. ([S])-®, we deduce the following proposition from Eqs. (jS)). 
Proposition 1. The triples I andl^ are given by 
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In any static, spherically symmetric space-time the geodesic equations imply that the vector L defined 
as L = —X X Hs a constant of the motion. The null geodesies considered here are assumed to be unbound. 
Consequently the magnitude of L is such that \L\ = li'm\x\^oo \^ x dx/cdt\ since i — > —(dx/cdt)oo when 
I a; I — > oo. So the quantity b defined by 

6=|-a;xT| (11) 

is the Euclidean distance between the asymptote to the ray and the line parallel to this asymptote passing 
through the center O as measured by an inertial observer at rest at infinity. Hence b may be considered 
as the impact parameter of the ray (see, e.g., Chandrasekhar 1983). Besides its geometric meaning, b 
presents the interest to be intrinsic, since it corresponds to a quantity which could be really measured. 



Substituting for 1^ from Eq. ([T0| into Eq. ([TT]) . introducing the zeroth-order distance of closest 
approach Vc defined as 

rc^-r—^ — lln^xn^l, (12) 
and then using (r^ + rs)\nA x Ugl/lxg - x^l = \N,is x + \Nab x tIbI, we get 

(7 + l)m I AT^B x I + \Nab x | 
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Using this expansion of &, we obtain the proposition which follows. 

Proposition 2. In terms of the impact parameter h, the triples and 1^ may be written as 
m\NAB X n,i| 
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3. DEFLECTION OF A LIGHT RAY EMITTED AT INFINITY 

Assume now that the ray arriving at Xb is emitted at infinity in a direction defined by a unit vector 
Nf.. Substituting N^. for Nab and —N^. for Ua in Eq. (|15p yields the expression of ^b, where 6 is furnished 
by the limit of Eqs. (IT^ and when —> cx) and — > — Afg. We can set a proposition as follows. 



Proposition 3. For a light ray emitted at infinity in a direction Ng and arriving IS given 



by 
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where PB{Nf,) is the unit vector orthogonal to defined as 
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and b is the impact parameter of the ray, namely 
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with rr = ralNf, x rial. 



The deflection of the ray at point Xb niay be characterized by the angle Axb made by the vector Ne 
and a vector tangent to the ray at Xb- We have 
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Substituting for 1^ from Eq. (IT5|) into Eq. (|19l) , and then introducing the angle (j)B between TVe and 
Ub defined by 

iVe-ne = COS0B, O<0B<7r, (20) 



we get 
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where the impact parameter given by Eq. (jl8l) may be rewritten as 
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It may be seen from the formulas given in Teyssandier & Le Poncin-Lafitte 2006 that ips + ^Xb is the 
angular distance between the center O and the source at infinity as measured at Xb by a static observer, 
i.e. an observer at rest with respect to the coordinates x^. It will be shown in a subsequent paper that 
this property implies that A^b can be regarded as an intrinsic quantity. 

The l/c^ term in Eq. (|21l) is currently used in VLBI astrometry. If b is replaced by its coordinate 
expression (|22p . it may be seen that Axb is given by an expression as follows 
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For a ray grazing a mass M of radius rg, the underbraced term in the r.h.s. of Eq. ((23|) generates 
a post-post-Newtonian contribution [Ax^g] grazing ~ —4(7-1- l)2(GAf/c2ro)^(rB/ro) which can be great 
if Tb ^ tq. For Jupiter, [Ax^b) grazing — 16.1 /ias if the observer is located at a distance from Jupiter 

= 6 AU: this value is appreciably greater than the level of accuracy expected for Gaia. However, this 
'enhanced' term is due to the use of the coordinate-dependent quantity instead of the intrinsic impact 
parameter b. This result confirms the conclusion recently drawn in Klioner & Zschocke 2010. 



4. CONCLUSION 

Deriving the second-order terms in the propagation direction of light from the time transfer function 
rather than from the null geodesic equations is a very elegant and powerful procedure. The application 
of this method to a ray emitted at infinity and received by a static observer located at a finite distance 
from the central mass is easy and yields an intrinsic characterization of the gravitational bending of light. 
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